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Abstract 

A measurement of the proton structure function F2(x,Q 2 ) is presented in the kinematic 
range 0.045 GeV 2 < Q 2 < 0.65 GeV 2 and 6- 10~ 7 < x < 1- 10~ 3 . The results were obtained 
using a data sample corresponding to an integrated luminosity of 3.9 pb _1 in e + p reactions 
recorded with the ZEUS detector at HERA. Information from a silicon-strip tracking 
detector, installed in front of the small electromagnetic calorimeter used to measure the 
energy of the final-state positron at small scattering angles, together with an enhanced 
simulation of the hadronic final state, has permitted the extension of the kinematic range 
beyond that of previous measurements. The uncertainties in F2 are typically less than 
4%. At the low Q 2 values of the present measurement, the rise of F2 at low x is slower 
than observed in HERA data at higher Q 2 and can be described by Regge theory with a 
constant logarithmic slope d\a.Fi/d\a.{l/x). The dependence of F2 on Q 2 is stronger than 
at higher Q 2 values, approaching, at the lowest Q 2 values of this measurement, a region 
where F2 becomes nearly proportional to Q 2 . 
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1 Introduction 



A remarkable feature of the proton structure function F2(x,Q 2 ) is its rapid rise at low x, observed 
by the HI and ZEUS collaborations [jjj at HERA. First observed for Q 2 values above lOGeV 2 , the 
persistence of this rise down to small Q 2 |2|-Q challenges our understanding of QCD. In a recent 
publication jl, the ZEUS collaboration has discussed the transition from deep inelastic scattering to 
photoproduction. It was found that standard non-perturbative approaches which apply in photopro- 
duction fail to describe the data in the region above Q 2 = 0.9 GeV 2 . Next-to-leading-order QCD fits 
are successful when taken down to these low Q 2 values. As Q 2 approaches 1 GeV 2 , however, these fits 
yield vanishing gluon densities at low x, while the sea-quark density remains finite. Such a "valence- 
like" gluon distribution, vanishing as x — > 0, seems unnatural even at low Q 2 and has led to much 
discussion ||. Precise measurements at low Q 2 are important in elucidating this subject. 

This letter presents a measurement of F 2 at low Q 2 (0.045 GeV 2 < Q 2 < 0.65 GeV 2 ) and low x 
(6 • 10~ 7 < x < 1 • 10~ 3 ). The data used correspond to an integrated luminosity of 3.9 pb -1 and were 
taken with dedicated triggers during six weeks of e + p running in 1997. Compared to a previous result^ 
[||, the new measurement covers a larger kinematic region with an improved statistical precision and 
systematic accuracy. This was made possible by the addition of a Beam Pipe Tracker in front of the 
Beam Pipe Calorimeter used for measuring the energy of the final-state positron at small scattering 
angles, and by an enhanced simulation of the hadronic final state. 



2 Kinematic variables and cross sections 

Inclusive deep inelastic positron-proton scattering, e + p — ► e + X, can be described in terms of two 
kinematic variables, x and Q 2 , where x is the Bjorken scaling variable and Q 2 the negative of the square 
of the four-momentum transfer. They are defined as Q 2 = —q 2 = — (k — k') 2 and x = Q 2 /(2P ■ q), 
where k and P are the four-momenta of the incoming positron and proton, respectively, and k' is the 
four-momentum of the scattered positron. The fractional energy transferred to the proton in its rest 
frame, y, is related to x and Q 2 by Q 2 = sxy, where s = 4E e E p is the square of the positron-proton 
center-of-mass energy. Here, E e = 27.5 GeV and E p = 820 GeV are the positron and proton beam 
energies, respectively. 

At the low Q 2 values of this measurement, where the contribution from Z exchange is negligible, 
the double-differential cross section for inelastic e + p scattering can be written in terms of F2 and the 
longitudinal structure function Fl as 

d 2 a lira 2 /, , , 0s \ , 

s^-^(< 1 + < 1 -»>*>«-» > ft)<i+«. w 

The QED radiative correction, S r , is a function of x and Q 2 , but, to a good approximation, independent 
of both F 2 and F L f§. 



3 Experimental setup and kinematic reconstruction 

The ZEUS detector has been described in detail previously In the present analysis, the scattered 
positron was detected in the Beam Pipe Calorimeter (BPC) and Beam Pipe Tracker (BPT) JTl], [12]j. 
The BPC was installed in 1995 to enhance the acceptance of the ZEUS detector for low-Q 2 events, 
where the positron is scattered through a small angle, and was used for a previous measurement of 
i<2 0. I n 1997, the BPT was installed in front of the BPC to complement the calorimetric energy 
measurement with precise tracking information. 

1 When the previous measurement of F2 at low Q 2 was published, it was found that existing models and parameteri- 
zations || were unable to describe the data. Subsequently, there have been various new approaches, as well as updates 
and improvements of existing ones, to describe -£2(2, Q 2 ) in the entire Q 2 range Q,^). The validity of these recent models 
is not discussed here. 
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3.1 Beam Pipe Calorimeter and Beam Pipe Tracker 

The BPC is a tungsten-scintillator sampling calorimeter with the front face located at Z = —293.7 cm, 
the center at Y = 0.0 cm, and the inner edge of the active area at X = 4.4 cm, as close as possible 
to the rear beam pipeQ. The BPC has an active area of 12.0 x 12.8 cm 2 in X x Y and its depth in 
Z corresponds to 24 Xq. The relative energy resolution as determined in test-beam measurements 
with 1-6 GeV electrons is AE/E = 11%/ y/E [GeV]. The scintillator layers are alternately subdivided 
into vertical and horizontal strips, 7.9 mm in width, referred to as X and Y fingers, respectively. The 
stacks of fingers with common X and Y positions are read out together by a wavelength shifter bar. 

The BPT is a silicon microstrip tracking device. In 1997, it was equipped with two detector planes 
to measure the X coordinate. The planes, of dimensions 6 x 6 x 0.03 cm 3 , have 576 strips with a 100 fim 



pitch. The detectors are read out with binary electronics 13]. The number of readout channels per 
detector is 512; 128 strips at the outer edge of a plane are read out in pairs, thus giving an effective pitch 
of 200 /im. The detector planes and front-end electronics are supported in a carbon-fiber composite 
structure attached to the front face of the BPC. The planes are positioned at Z = —252.7cm and 
Z = —279.1 cm and have centers at X = 6.85 cm and Y = 0.0 cm. 

Scattered positrons traveling from the interaction point toward the BPC and BPT leave the beam 
pipe through a 1.5 mm thick (1.6% Xq) aluminum exit window positioned at Z = —249.8cm. The 
fiducial area of the measurement was defined by the overlap of the exit window with the acceptance 
of BPC and BPT, and was thus limited to a D-shaped region in the range 5.2 cm < X < 9.3 cm and 
—2.3 cm < Y < 2.8 cm at Z = —293.7 cm. Only positrons scattered into the fiducial area were used in 
the analysis, yielding an angular acceptance of 18mrad < 6 e < 32mrad for particles emanating from 
the nominal interaction point. 

3.2 Measurement of the scattered positron 

The combination of calorimetric and tracking information for the measurement of the four-momentum 
of the scattered positron resulted in significant improvements with respect to the previous analysis. 
The BPT tracking information proved crucial in suppressing beam- and neutral-particle-related back- 
grounds. It also allowed a better control of the systematic uncertainties in the position-dependent 
corrections to the BPC energy measurement, the development of an improved position reconstruc- 
tion algorithm for the BPC, a direct determination of the fiducial area of the measurement, and a 
significant reduction of the uncertainty in the position of the BPC. Moreover, the fact that the lon- 
gitudinal coordinate of the event vertex was reconstructed from the scattered-positron track reduced 
the dependence of the reconstruction of kinematic variables on the hadronic final state. 

The energy deposited in the BPC by the scattered positron was reconstructed as the sum over 
all energies measured in the BPC fingers. Calibration constants for individual stacks of fingers were 
determined using kinematic-peak events) 3 ], taking into account position-dependent corrections for trans- 
verse shower leakage from the BPC, light attenuation in the scintillator fingers, and the non-uniformity 



caused by the gaps between fingers [12]. After the calibration procedure, the BPC energy response 
was found to be uniform to within 0.3%, and the absolute energy scale at 27.5 GeV was known to 
0.3%. The linearity of the energy response was estimated through a simulation based on the results 
of a scan with a 60 Co source; it was found that the energy scale at 4 GeV was accurate to within 1% 



o 



The results of the simulation were confirmed by a study of QED Compton scattering events [|l^] . 
The transverse position of a positron shower in the BPC was reconstructed by using the correlation 
between the shower position relative to the center of the struck stack of BPC fingers and the energy 



2 The ZEUS right-handed Cartesian coordinate system has its origin at the nominal interaction point, the Z axis 
pointing in the proton beam direction (referred to as forward direction), and the X axis pointing toward the center of 
HERA. The polar angle ■& is measured with respect to the positive Z axis; it is convenient to also define 8 — n — The 
pseudorapidity is defined as r\ — — ln(tan(i9/2)). 



3, 



A cut yjB < 0.04 (see Section 3.3) selects events for which the energy of the scattered positron sharply peaks within 



2% of the beam energy, providing a good calibration source (kinematic-peak events). 



2 



deposited in its two neighbors. The function correlating the distribution of energy to the shower 
position was determined by comparing positron positions reconstructed separately with the BPC 



and the BPT, respectively [12|. The resulting average BPC position resolution was AX = Ay = 
0.22cm/V-E [GeVJ. 

The projection of the positron scattering angle on the X-Z plane, 9x, the positron impact point 
on the BPC front face in X, and the longitudinal position of the event vertex, Z vtx , were determined 
by reconstructing a BPT track as the straight line joining the hits in the two planes and assuming 
the average value for the X coordinate of the event vertex, determined with the Central Tracking 
Detector (CTD, see below). A study of simulated events showed that the effect of the magnetic field 
was negligible. The value of 9y could only be reconstructed from the shower position in the BPC. 
However, since for positrons scattered into the fiducial area 9x was significantly larger than 9y, the 



impact of the limited 9y resolution on the reconstruction of 9 e = y 9 X + 9 Y was small. On average, 
the angular resolution was A9 e = 0.2mrad, and the vertex resolution was A^ v ^x — 3 cm. 

The X positions of the BPC and BPT were determined with a minimization procedure based on 
the comparison between the event vertex reconstructed from the BPC and BPT information and that 



reconstructed from the CTD tracking information |12|. Several studies were performed to determine 



the stability of the result: the largest contribution to the systematic error comes from the uncertainty 
in the relative Z position of the BPT planes, estimated from a survey to be 300 /im. The resulting 
total accuracy in the X positions of the BPC and BPT is 200 /um. The tracking efficiency of the BPT 
was 90.4 ± 1.5%, the inefficiency being primarily due to dead strips. 

3.3 Measurement of the hadronic final state 

The hadronic final state was measured using the uranium-scintillator calorimeter (CAL) and the 
CTD. The CAL covers 99.7% of the total solid angle. Under test-beam conditions, its relative energy 
resolution is AE/E = 18%/y/E [GeVj for electromagnetic showers and AE/E = 35%/V-E [GeVJ 
for hadronic showers. The CTD operates in a 1.43 T solenoidal magnetic field, and has a relative 
resolution for full-length tracks of Apx/pr = 0.0058 • px © 0.0065 © 0.0014/p^, with px measured 
in GeV (the symbol © denotes addition in quadrature). The interaction vertex is measured with a 
typical resolution along (transverse to) the beam direction of 4 mm (1mm). 

Global variables characterizing the hadronic final state were measured by summing all recon- 
structed hadronic final-state objects. A combination of clusters of energy deposits in the CAL and the 



corresponding tracks measured in the CTD [16| was found to provide the best resolution in determining 
the quantities 



pf d = J2 E i Sin #i COS & » PY^ = E E i Sil1 * Sin & > =Y, E i COS * > 



E had =J2 E i, p T &d = V \px A f ' + (PY d Y \ ^ = Ehad - Pz^, 



^had 

Vm = 2E e 



(2) 



3.4 Kinematic reconstruction 

Two methods have been used for the reconstruction of the kinematic variables Q 2 and y; the variable 
x was then determined from x = Q 2 /(sy). The first method uses the energy, E' e , and angle, 9 e , of the 
scattered positron ("electron method") to determine 

Ql = 2E e E , e {l + cos# e ) « E e E' e 9 2 e , 



F' F' 

1-— f (l-costf e ) w 1--^. 

r) Tp \ c / TP 

L£j e tj e 



(3) 



The second method is used for low values of y, where the y resolution of the electron method can 
be improved by combining the scattered-positron variables with hadronic final-state variables ("eS 
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method' 



|l7| ) to determine 




Ql 



2E e 5 had 



(4) 



VeT. = 



5 ' 5 



where 5 is defined as 5 = <5 had + E' e (l — cos$ e ). Four-momentum conservation requires 5 = 2E e for 
fully contained and perfectly measured events. 

In the analysis, the kinematic variables of the events were reconstructed with the electron method 
for y e > 0.08, and with the eS method for y e -£ < 0.08. The resulting relative resolution in y ranged 
from about 5% at high y to 25% at low y. The relative Q 2 resolution was about 10% in the entire 
kinematic range of the present measurement. The use of the eS method allowed the extension of the 
measured kinematic range to y values as low as 0.005. 

3.5 Luminosity measurement 

The luminosity was measured from the rate of photons from the Bethe-Heitler process, ep — > ep'y, 
using a lead-scintillator calorimeter |l8| positioned at Z = —107 m to detect photons scattered through 
angles of less than 0.5 mrad. The data set used in the analysis corresponded to an integrated luminosity 
of 3.9 pb -1 , with an uncertainty of 1.8%. 

4 Trigger, event selection and background 

The event selection was based mainly on the requirement of a well-reconstructed positron in the BPC 
and in the BPT, while additional cuts on the hadronic final state suppressed background and limited 
effects of resolution smearing (event migrations) and radiative corrections. 

Events were selected online by the ZEUS three-level trigger system. The trigger required a mini- 
mum energy deposit in the BPC, a timing compatible with an ep interaction, and imposed requirements 
on energy deposits from the hadronic final state in the CAL. No BPT information was used in the 
trigger. 

The offline event selection imposed tighter requirements than those applied in the trigger. Using 
an independently triggered control sample, the trigger was found to be more than 99.5% efficient for 
events passing the offline cuts. 

The scattered positron was required to have an energy, E' e , of at least 4.4 GeV for events recon- 
structed with the electron method, or at least 20 GeV for the eS method. The positron position at 
the BPC front face as extrapolated from the BPT measurement had to lie within the fiducial area. 
In order to identify electromagnetic showers and reject hadrons, the transverse size (energy-weighted 
r.m.s.) of the shower in the BPC was required to be less than 0.8 cm. In order to limit migrations from 
low y, a cut yjB > 0.06 was applied for events reconstructed with the electron method (yjB > 0.004 
for the e£ method). 

Background came primarily from events in which a particle from the hadronic final state, usually 
a photon from the decay of a 7r°, produced an electromagnetic shower in the BPC, while the scattered 
positron escaped undetected. Such events originated mostly from photoproduction (with Q 2 lower 
than measured in the present analysis). In addition, photons from initial- and final-state radiation 
(ISR and FSR) could produce an electromagnetic shower in the BPC. Since the longitudinal mo- 
mentum carried away by the undetected particles (usually escaping through the rear beam-pipe hole) 
decreases 5, a cut 5 > 30 GeV was applied to reject both photoproduction background and radiative 
events with a hard initial-state photon. In addition, background from neutral particles was suppressed 
through the requirement that the extrapolated track from the BPT match the position of the BPC 
shower in X within five times the average BPC position resolution. The remaining contamination was 
determined from the number of events in which the scattered positron was detected in an electromag- 
netic calorimeter positioned at Z = —35 m, designed to tag positrons with Q 2 < 0.01 GeV 2 scattered 
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through very low angles, scaled by the calorimeter tagging efficiency (estimated using a sample of sim- 
ulated photoproduction events passing all selection cuts). The background contamination was found 
to be less than 1.5% at high y and negligible at low y, and was statistically subtracted from the data 
sample. The uncertainty associated with this background correction was estimated by comparing it 
to the number of simulated photoproduction events (in a luminosity-normalized sample) which passed 
all analysis cuts and taking the difference to be the uncertainty. 

In order to suppress beam-related background, the longitudinal coordinate of the event vertex as 
reconstructed by the BPT was required to lie within 90 cm of the nominal interaction point. A cut 
5 < 65 GeV rejected events from processes other than ep interactions. The fraction of such background 
events, determined from events recorded at times when there was no ep beam crossing, was found to 
be below 0.5% and was neglected. 



5 Analysis 

5.1 Monte Carlo simulation 

Monte Carlo (MC) simulations were used to characterize the accuracy of the kinematic reconstruction, 
to determine the efficiency of selecting events, to estimate the background rate, and to extract i*2- 
Non-diffractive processes including first-order QED radiative corrections were simulated using the 
HERACLES 4.6.1 program with the DJANGOH 1.1 interface 111 to the QCD programs. The program 



RAPGAP 2.06 [20] was used to simulate diffractive processes in which the incoming proton emits a 
Pomeron that has partonic structure. The color-dipole model of ARIADNE 4.08 |2lj] was used for 
the simulation of the parton shower, while the hadronization was based on the Lund string model of 
JETSET 7.410 §§. Photoproduction background events were simulated using PYTHIA 5.724 @. 
All generated events were passed through a detector simulation based on GEANT 3.13 [p4|]. 

In order to improve the simulation of the hadronic final state, the diffractive and non-diffractive 
samples were mixed in a proportion determined from the data (while maintaining a fixed overall 
MC normalization). This was crucial at high y and low Q 2 , where the trigger and offline event 
selection efficiencies for diffractive and non-diffractive events were significantly different, owing to the 
different topologies of the hadronic final state. The mixing ratio was determined by optimizing the 
agreement between data and MC in distributions of hadronic final-state quantities according to a 
least-squares minimization procedure applied to histograms of these quantities fi~2|| . The hadronic 
variables chosen for the determination of the mixing ratio, parameterized^ as a function of x, were the 
pseudorapidity of the most forward hadronic energy deposit or track, 77 max , and the variables <5 had and 
Pj? ,d in combination. The average of the two results for the mixing fraction was used for the extraction 
of i*2, while the difference between them was taken into account as a contribution to the systematic 



uncertainties (see Section 5.4). The fraction of diffractive events used for the extraction of F2 varied 
between 10% and 25% over the x range of the measurement. 

Figure [l] shows comparisons between distributions of reconstructed positron and hadronic final- 
state variables in data and simulation for 0.06 < yjB and y e < 0.74. Figure |^a-c shows comparisons 
between distributions of reconstructed kinematic variables. Good agreement between data and simu- 
lation was obtained in this region and the agreement is of similar quality in the remaining kinematic 
region used for the extraction of F<i- 

5.2 Binning of the data 

Figure §d shows the bin boundary grid used to extract F%. The bin widths in y are approximately 
three times the y resolution at low y, and twice the resolution at higher y. The bin widths in Q 2 are 

4 Although in principle the mixing ratio between the diffractive and non-diffractive Monte Carlo samples may depend 
on both x and Q 2 , it was found that it could be parameterized with sufficient accuracy as a function of x only. 
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about four times the Q 2 resolution. In the region of overlap, the bin boundaries were chosen to be 
identical to those used in the previous measurement [Q. 

The geometrical acceptance (defined as the fraction of events where the scattered positron hits the 
fiducial area) was determined from the Monte Carlo simulation to be between 4% and 14% for most 
of the bins, decreasing to 1% for a few bins. The acceptance (defined as the fraction of events passing 
all selection cuts and generated in a bin that were also reconstructed in that bin) varied between 50% 
and 70%, and the purity (defined as the fraction of events passing all selection cuts and reconstructed 
in a bin that were also generated in that bin) typically ranged from 40% to 65%, decreasing to 30% 
at the lowest y. 



5.3 Extraction of F 2 

The proton structure function F 2 was measured using an iterative bin-by-bin unfolding method. The 
measured value of F 2 at a given point (x,Q 2 ) inside a bin was obtained as the ratio of the numbers 
of events reconstructed in the bin in data and MC, multiplied by the corresponding F 2 value of the 
parameterization used in the MC sample. The point (x,Q 2 ) at which F 2 is quoted was chosen such 
that y and Q 2 are round numbers close to the center of the bin. 

First-order QED radiative corrections were included in the MC simulation used for the unfolding. 
The uncertainty in these corrections is due principally to the uncertainty in F 2 at x and Q 2 values 
outside the region of this measurement, which was estimated to be 6% from an analysis of events with 



ISR 1 12] and was taken into account in the systematic error calculation. Higher-order QED radiative 



corrections, including soft-photon exponentiation, were evaluated using the program HECTOR 25] in 
the leading-log approximation. They were found to be less than 0.2% and were thus neglected. 

Since the bin- by-bin unfolding requires a precise simulation of migration effects, which in turn 
requires F 2 to be approximately the same in data and MC, the MC sample was iteratively re-weighted 
by the extracted F 2 . The contribution from Fl was neglected in this procedure. The MC events 
were weighted according to the ALLM97 parameterization in the first iteration, and according to 
a function of the forrn^] 




in further steps. At each iteration, oir was fixed to the value 0.5 and the values of the remaining four 
parameters in Eq. (0) were obtained from a fit to the results obtained in the previous step; in order to 
constrain the fit at high x, photoproduction data [28 were also included^. The F 2 results converged 



after a few iterations. 

To correct for the effect of Fl, the value of R = Fl/(F 2 — Fl) from the BKS model [^3] was 
used, which, in the range of interest here, can be parameterized to a very good approximation by 
R = 0.165 • Q 2 /m 2 p , where m p = 0.77 GeV is the p meson mass. This changed the extracted F 2 in 
the bins of the present measurement by at most 3% with respect to the values determined assuming 
F L = 0. 



5.4 Systematic uncertainties 

The systematic uncertainties in the measured F 2 values were determined by studying the stability 
of the results under variations of the reconstruction parameters entering detector calibration and 

5 This phenomenological parameterization is based on the combination of a simplified version of the generalized vector 
meson dominance model j26| for the description of the Q 2 dependence and Regge theory j27) for the description of the 
x dependence of Fi. 

6 For this purpose, Eq. ^ was re-written using a^ ot p = (An 2 a/Q 2 )F2 and W 2 = Q 2 /x, where W denotes the photon- 
proton center-of-mass energy. 
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alignment, of the simulated detector efficiency, of the event selection cuts, and of several other aspects 
of the analysis procedure. The reconstruction parameters were varied within their uncertainties. The 
variations of the selection cuts took into account the degree of arbitrariness in the cut values. The 
following checks were performed (the effect on F2 is given in parentheses): 

• Analysis cuts: 

— variation of the cut 5 > 30 GeV by ± 2 GeV (mostly below 1.5%, up to 4% at high y); 

— variation of the cut y JB > 0.06 by ±0.01, and of the cut y JB > 0.004 by ±0.001 (up to 1.5% 
at medium y); 

— variation of the BPC shower-width cut at 0.8 cm by ±0.1 cm (up to 1.5% at low positron 
energies) ; 

— variation of the 5cr BPC shower/BPT track-match cut by ±2<t (up to 1.5% at low positron 
energies) ; 

— change of the BPT vertex cut from ±90 cm to ±50 cm (up to 3.5% in bins with low geo- 
metrical acceptance); 

— variation of the fiducial cut in X by ±1 mm (mostly below 1.5%, up to 6% in bins with low 
geometrical acceptance); 

— variation of the fiducial cut in Y by 11mm (up to 2.5%). 

• Detector calibration, alignment, and efficiency: 

— variation of the BPC energy scale by ±0.3% (up to 2.5% at low y); 

— variation of the BPC energy response by ±1% at 4 GeV, decreasing linearly to 0% at 
27.5 GeV (up to 2.5% at medium and high y); 

— variation of the absolute BPC/BPT position in X by ±200 /xm (mostly 1.5%, up to 3.5% 
in some bins); 

— variation of the simulated BPT efficiency by ±1.5% (1.5%); 

— variation of the CAL energy scale for hadrons by ±3% (up to 2.5% at high and low y). 

• Analysis procedure: 

— variation of the fractions of DJANGO and RAPGAP events in the MC sample by half of 
the difference between fractions obtained from separate optimizations in 77 max or in 8 had 
and plj? d (mostly below 1%, up to 8% in the two highest y bins); 

— variation of the subtracted photoproduction background by ±200% and —100% (up to 2.5% 
at high y); 

— re-weighting of F% in the MC sample, outside the bins used for the measurement, by ±6%, 
which induces a variation of the number of events migrating into the measurement region^] 
(up to 3% in some bins). 

The total systematic error was computed as the quadratic sum of the individual contributions, 
separately for positive and negative deviations. The average statistical error is 2.6% and the average 
systematic error 3.3%. In most bins, the systematic error has a magnitude similar to that of the 
statistical error. No individual contribution to the systematic error dominates, except in the two 
highest y bins, where the error is dominated by the uncertainty in the fraction of diffractive events. 
An additional overall normalization uncertainty of 1.8% due to the luminosity measurement is not 
included in the systematic error. 

7 This check is particularly sensitive to the magnitude of radiative corrections due to migrations of unrecognized ISR 
events from higher x and lower Q 2 , but also to the amount of migrations due to detector resolution effects. 
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6 Results 



The measurement of F% presented here uses data in the kinematic region 0.04 GeV 2 < Q 2 < 0.74 GeV 2 
and 5.3 • 10~ 7 < x < 1.6 • 10~ 3 , corresponding to 0.005 < y < 0.84. The values of F2 extracted in 70 
bins are listed in Table Fi gure |3j shows these F2 values as a function of x for different bins of Q , 
together with previous ZEUS and HI data at low Q 2 and with data at higher x from the fixed- 

target experiment E665 |}0|. The curve denoted as "ZEUS Regge fit" represents the parameterization 
of Eq. (H|), with the values of the parameters resulting from the last iteration in the extraction of F^: 
A u = 147.8 ± 4.6 fib, a n = 0.5, A w = 62.0 ± 2.3/xb, a P = 1.102 ± 0.007, M 2 = 0.52 ± 0.04 GeV 2 
(statistical and systematic errors added in quadrature). These results are in good agreement with 
those previously obtained 0. 

The errors in F2 are significantly reduced compared to the previous measurement Q. While in 
general good agreement is observed in the region of overlap, the present results are slightly lower than 
the previous ones at the lowest Q 2 and low x. At the highest x of this measurement, the present 
data overlap with data from E665. Although the F2 values do not differ significantly in the region of 
overlap, if the ZEUS Regge fit is extrapolated to higher x into the E665 region, it is found to lie about 
15% above the E665 values. 

The rise of the proton structure function F2 at low x, observed to be steep in HERA data at 
higher Q 2 , persists down to small Q 2 , but becomes shallower as Q 2 decreases into the range of the 
present measurement. At the low Q 2 values of this measurement, the rise of F2 at low x is well 
described by Regge theory assuming a constant logarithmic slope dlnFz/dln^l/x), as reflected in the 
good agreement between the data and the ZEUS Regge fit. 

Figure |] shows F2 as a function of Q 2 for different bins of y, together with previous ZEUS and 
HI data |],|,§. At higher Q 2 , F 2 is roughly independent of Q 2 . The curve denoted as "ZEUS QCD 
fit" |H illustrates that this behavior is well described by next-to-leading-order QCD fits down to Q 2 
values of about 1 GeV 2 . On the other hand, it is clear that F2 must acquire a stronger Q 2 -dependence 
at sufficiently low Q 2 , since conservation of the electromagnetic current requires F2 to vanish like Q 2 
as Q 2 — > 0. Dynamical mechanisms, e.g. parton saturation at small x values [3l| , can produce such a 
behavior at low Q 2 . The present data exhibit a smooth transition to a stronger Q 2 dependence in the 
Q 2 range between 0.1 GeV 2 and lGeV 2 , approaching, at the lowest Q 2 values of this measurement, a 
region where F2 becomes nearly proportional to Q 2 . This transition is well described by the ZEUS 
Regge fit. 

7 Summary 

The proton structure function F2(x,Q 2 ) has been measured in the kinematic range 0.045 GeV 2 < 
Q 2 < 0.65 GeV 2 and 6 • 10~ 7 < x < 1 • 10~ 3 , using an e + p data sample corresponding to an integrated 
luminosity of 3.9 pb -1 . The addition of a Beam Pipe Tracker in front of the Beam Pipe Calorimeter 
and an enhanced simulation of the hadronic final state have resulted in coverage of a larger kinematic 
region and in improved statistical precision and systematic accuracy compared to previous results 
obtained with the Beam Pipe Calorimeter alone. 

At the low Q 2 values of the present measurement, the proton structure function F2 rises more 
slowly with x than observed in HERA data at higher Q 2 . This slow rise can be described by Regge 
theory with a constant logarithmic slope d\nF2 / d\n(l / x) . Furthermore, the F2 data presented here 
exhibit a stronger Q 2 dependence than observed at higher Q 2 , approaching, at the lowest Q 2 values 
of this measurement, a region where F2 becomes nearly proportional to Q 2 . 
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ZEUS 1997 



• ZEUS BPT 1997 □ MC DJANGO □ RAPGAP 




Figure 1: Distributions of reconstructed quantities in data and MC simulation in the region denned by 
0.06 < ?/jb and y e < 0.74: a) energy, E' e , of the positron measured in the BPC; b) scattering angle, 6 e , 
of the positron measured in the BPT; c) p^ ad ; d) (5 had ; e) r/ max ; f) 5. The points denote measured data, 
the light shaded histogram is the sum of non-diffractive (DJANGO) and diffractive (RAPGAP) MC, 
the dark shaded histogram and the dashed line represent the individual contributions from RAPGAP 
and DJANGO, respectively. 
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Figure 2: a)-c) Distributions of reconstructed kinematic quantities in data and MC simulation in the 
region defined by 0.06 < yje and y e < 0.74: a) Q 2 ; b) x e ; c) yje- The points denote measured data, 
the light shaded histogram is the sum of non-diffractive (DJANGO) and diffractive (RAPGAP) MC, 
the dark shaded histogram and the dashed line represent the individual contributions from RAPGAP 
and DJANGO, respectively, d) The bins in the kinematic plane (Q 2 vs. x). The bin boundaries in 
y are 0.005, 0.01, 0.02, 0.04, 0.08, 0.16, 0.23, 0.30, 0.37, 0.45, 0.54, 0.64, 0.74, 0.84; those in Q 2 are 
0.040, 0.055, 0.075, 0.10, 0.13, 0.17, 0.21, 0.27, 0.35, 0.45, 0.58, 0.74 GeV 2 . The border at y = 0.08, 
above which the electron method and below which the eS method was used to reconstruct the event 
kinematics, is also indicated. 
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Figure 3: Measured F2 vs. x in bins of Q 2 . The data from the present measurement are indicated by 
filled circles. The solid line shows the ZEUS Regge fit. Open circles denote the results from a previous 
analysis, filled and open triangles denote other measurements from ZEUS and HI, respectively, and 
squares denote results from E665. These other measurements have been shifted to the Q 2 values of 
the present measurement using the ALLM97 parameterization. The inner error bars represent statis- 
tical errors, the outer ones the sum in quadrature of statistical and systematic errors; normalization 
uncertainties are not included. 
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Figure 4: Measured F2 vs. Q 2 in bins of y. The data from the present measurement are indicated 
by filled circles. Triangles, filled squares and open circles denote other measurements from ZEUS and 
HI. The data have been scaled by the numbers in parentheses for clarity of presentation. The solid 
line at low Q 2 shows the ZEUS Regge fit, the dashed line at higher Q 2 the ZEUS QCD fit. The 
other measurements have been shifted to the y values of the present measurement using the ALLM97 
parameterization. 
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